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Monte Carlo (MC) analysis of the Goldstone mode singularities for the transverse and the longitudinal correla- 
tion functions, behaving as G±[k) - ak~^ L and G||(k) = bk~^t in the ordered phase at k — > 0, is performed 
in the three-dimensional O(n) models with n = 2,4, 10. Our aim is to test some challenging theoretical pre- 
dictions, according to which the exponents Aj_ and Ay are non-trivial (3/2 < Aj_ < 2 and < X\\ < 1 in three 
dimensions) and the ratio bM 2 1 a 2 (where M is a spontaneous magnetization) is universal. The trivial standard- 
theoretical values are Aj_ = 2 and A|| =1. Our earlier MC analysis gives Aj_ = 1.955 ± 0.020 and Ay about 0.9 
for the 0(4) model. A recent MC estimation of An , assuming corrections to scaling of the standard theory, yields 
A|| = 0.69±0.10 for the 0(2) model. Currently, we have performed a similar MC estimation for the 0(10) model, 
yielding Aj_ = 1.9723(90). We have observed that the plot of the effective transverse exponent for the 0(4) 
model is systematically shifted down with respect to the same plot for the 0(10) model by AAx = 0.0121(52). 
It is consistent with the idea that 2 - Aj_ decreases for large n and tends to zero at n — oo. We have also veri- 
fied and confirmed the expected universality of bM 2 la 2 for the 0(4) model, where simulations at two different 
temperatures (couplings) have been performed. 
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1. Introduction 

Our work is devoted to the Monte Carlo (MC) investigation of the Goldstone mode effects in n-com- 
ponent vector-spin models (0(n) models), which have 0{ri) global rotational symmetry at zero external 
field h. The Hamiltonian J€ is given by 



T F 



\<ij) 



(1.1) 



where T is temperature, s; = s(x,) is the n-component vector of unit length, i. e., the spin variable of the 
i-th lattice site with coordinate Xj, and p is the coupling constant. The summation takes place over all 
nearest neighbors in the lattice with periodic boundary conditions. 

The Fourier-transformed longitudinal and transverse correlation functions are 



G,|(k) = AT^GnWe- 

X 

G x (k) = AT^GiWe- 



ikx 



ikx 



(1.2) 
(1.3) 



where G\\ (x) and G± (x) are the corresponding two-point correlation functions in the coordinate space. 
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In the thermodynamic limit below the critical temperature (at p 1 > /3 C ), the magnetization M{h) and 
the correlation functions exhibit Goldstone mode power-law singularities: 

M{h) - M(+0) oc hP at h^O, (1.4) 
Gx(k) = afc _Al at fc = +0andfc^0, (1.5) 
G||(k) = fofc" A i at /i=+0andfc^0, (1.6) 

where a and b are the amplitudes. 

There exist different theoretical predictions for the values of the exponents in these expressions. In 
a series of theoretical works (e. g., (lrtZD), it has been claimed that these exponents are exactly p = 1/2 at 
d — 3, Aj_ = 2 and Ay = 4- d. Here, d is the spatial dimensionality 2 < d < 4. These theoretical approaches 
are further referred to as the standard theory. 

More non-trivial universal values are expected according to |8], such that 

rf/2<A_L<2, (1.7) 
A|| = 2A ± - <i , (1.8) 
p = (d/A x )-l (1.9) 

hold for 2 < d < 4. These relations were obtained in [8] by analyzing self-consistent diagram equations 
for correlation functions without cutting the perturbation series. As introduced in |9, 10], we will call 
this approach the GFD (grouping of Feynman diagrams) theory. Apart from the mathematical analysis, 
certain physical arguments were also provided |8] to show that X± — 2 could not be the correct result for 
the XY model (n = 2) within 2 < d < 4. 

Several MC simulations were performed in the past fllUlHl to verify the compatibility of MC data with 
some standard-theoretical expressions, where the exponents are fixed. In recent years, we performed a 
series of accurate MC simulations (ljj[lJi[l6j] for remarkably larger lattices than previously with an aim 
to reexamine the theoretical predictions by evaluating the exponents in <1.7> - fL9l . In particular, lattices 
of the linear sizes L =£ 512 for n — 2 and L =£ 350 for n — 4 were simulated in our papers (^llH and 1 16], 
respectively. These L values remarkably exceed the largest sizes simulated by other authors, i. e., L = 160 
for n - 2 in (l3|] and L = 120 for n = 4 in fl2Lfl3l . In the current work, the O(10) model is simulated up to 
L = 384. 

The relations dl.7> and il.8) are consistent with MC simulation results for the 3D 0(4) model 116], 
where an estimate Ax = 1.955 + 0.020 was found. It was also stated that the behavior of the longitudinal 
correlation function is well consistent with X\\ about 0.9 rather than with the standard-theoretical value 
A|| = 1. According to dl.9t . we have 1/2 < p < 1 in three dimensions. It is consistent with the MC estimate 
p = 0.555(17) for the 3D XY model (Hi, which corresponds to A ± = 1.929(21) according to O). A clear 
MC evidence that the behavior of G\\ (k) is not quite consistent with the standard-theoretical predictions 
has been recently provided 1 10], where an estimate X\\ = 0.69 + 0.10 has been obtained for the 3D XY 
(i. e., 3D 0(2)) model (at (5 = 0.55), assuming corrections to scaling of the standard theory. 

In the actual study, we have extended our MC simulations and analysis to include the n - 10 case 
and to test the n-dependence of the exponents. Apart from the exponents, we have performed here an 
extended analysis of the 0(4) model to verify the expected universality of the ratio bM 2 1 a 2 Q], where 
M = M(+0) is the spontaneous magnetization, a and b are the amplitudes in dl.St and dl.6> . 

2. Simulation results 

We simulated the 3D O(10) model by a modified Wolff cluster algorithm, used also in (l^[l^]. an d 
evaluated the Fourier-transformed correlation functions by techniques described in 1 16] . The standard 
Wolff cluster algorithm 1 17] was modified to enable simulations at nonzero external field h. Simple cubic 
lattices of the linear size up to L — 384 were simulated at /3 — 3 and h=\h\= h m i n , 2h m i n , 4h m [ n , where 
^min = 0.00021875. The coupling constant p 1 = 3 corresponds to the ordered phase, since the spontaneous 
magnetization M(+0) is about 0.467 in this case — see section[5]for details. This value of M(+0) is com- 
parable with those for the 0(2) and 0(4) models in our previous MC simulations (l5l[l6j]. The simulation 



43005-2 



Goldstone mode singularities in 0(«) models 




Figure 1. Log-log plots of the transverse correlation function Gj_(k) at /z = /z m i n = 0.00021875 and L = 350 
(solid circles), h = ft m ; n and L = 256 (pluses), h = 2fa m i n and L = 384 (empty circles), as well as at h = 
4/zjnin and L = 384 (empty diamonds). Statistical errors are about the symbol size or smaller. The lower 
value k* of the wave vector magnitude, used in estimations of the exponent X±, is indicated by a vertical 
dashed line. 



results for the correlation functions Gj_(k) and G||(k) in the (100) crystallographic direction at the three 
values of h and different sizes L are illustrated in ngures[l]and[2] It is important for an estimation of the 
exponents Aj_ and X\\ to ensure that the finite-size as well as finite-^ effects are small. This condition is 
satisfied for k> k* , where the values of k* are indicated in the figures by vertical dashed lines. 




Figure 2. Log-log plots of the longitudinal correlation function G||(k) at h = fa m j n = 0.00021875 and L = 
350 (solid circles), h = /z m j n and L = 256 (pluses), h = 2h m [ n and L = 384 (empty circles), as well as at 
h = 4fa m j n and L = 384 (empty diamonds). Statistical errors are about the symbol size. The lower value 
of the wave vector magnitude, used in estimations of the exponent X\\ , is indicated by a vertical dashed 
line. 
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3. Estimation of the exponents 

Here we estimate the exponents Aj_ and A \\ , describing the behavior of the correlation functions in the 
limit k — » 0, ft — ► 0, L — » oo, taking the limit L — «■ oo at first, followed by ft — » 0. For this purpose, first we find 
good approximations of the effective exponents at ft — » 0, L — ► oo, and then fit these fc-dependent effective 
exponents to evaluate their asymptotic values at k — ► 0. By comparing the simulation results for different 
L and ft, we conclude that the largest-L and smallest-ft data for k > k* with a good enough accuracy 
correspond to the thermodynamic limit at ft — +0, i. e., ft — * 0, L —> oo. We have tested this precisely by 
looking how the estimates of the effective exponents depend on L and ft. This method of analysis was 
applied in (l^ll^l- The effective transverse exponent A e ff(fc) for the 0(4) model was evaluated in flill 
from the slope of the InGx(k) vs lnfc plot within [k,4k]. Here we use a wider interval — [k,6k], because 
we have found that the A e ff(fc) data in this case can be perfectly fit by a parabola 

Aeff(fc) = Ax + a\k+ 0.2k 2 , (3.1) 

the finite-size and finite-ft effects being very small. The ansatz <3.1> is consistent with the general state- 
ment limfc^o A e ff (fc) = Ax (in the thermodynamic limit at ft = +0) and with corrections to scaling of the 
standard theory, where the correlation functions are supplied with correction factors in the form of an 
expansion in powers of k 4 ~ d and k d ~ 2 yj|5|]. Some of the fit results are shown in figure[3] We have per- 
formed a series of fits at different sizes L for the smallest-ft value ft = ft mm = 0.00021875. At the largest 
size L = 350 for this ft, the effective exponent A e ff(fc) was fit within k e [fcg, fes], where k( - 2n£l350 
are the possible discrete values of k. Similar fit intervals were chosen for all L. These fits to (3.1) give us 
A x = 1.9680(84) at L = 128, A ± = 1.9840(98) at L = 192, A ± = 1.9727(86) at L= 256 and A ± = 1.9723(90) 
at L = 350. As we can see, the finite-size effects are smaller than the statistical error bars. The A e ff(fc) 
data for L = 350 and L = 256 at ft = ft m i n are shown in figure [3] by solid circles and exes, respectively. 
The corresponding fit curves lie practically on top of each other. Therefore, only that one for L = 350 is 
shown by solid line. The fit curves for h = 2ft m j n and ft = 4ft m i n at L — 384 (the largest size) are also de- 
picted here to see the finite-ft effects. These fits give us X]_ = 1.9251(86) at ft = 4ft m i n and X± = 1.9666(91) 
at ft = 2ft m i n . A rather fast convergence to the ft = +0 limit is evident. According to this discussion, the 
fit result Ax = 1.9723(90), obtained at ft = ft m ; n and L = 350, with a good accuracy corresponds to the 
thermodynamic limit at h — +0. Besides, the systematical errors due to finite-size and finite-ft effects are 
probably smaller than the statistical error bars. 
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Figure 3. The transverse effective exponent k e ^(k), evaluated from the lnGx(k) vs In A: fits within [k,6k] 
at ft = ft m i n = 0.00021875 and L = 350 (solid circles), h = ft min and L = 256 (exes), ft = 2ft min and L = 384 
(empty circles), as well as at ft = 4ft m ; n and L = 384 (empty diamonds). The fits to (3. It of the largest-L 
data at ft = ft m j n , ft = 2ft mm and h = 4ft mm are shown by solid, dashed and dotted lines, respectively. 
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Another possible source of systematical errors is the existence of non-trivial corrections to scaling, 
which are not included in 13. H . These are corrections to scaling of the GFD theory |8], represented by 
an expansion in powers of k 2 ~ kL , k x t and Nevertheless, the actual estimation, where only the 
standard-theoretical corrections have been included, is well justified as a test of consistency of the stan- 
dard theory. The existence of a small correction-to-scaling exponent 2 - Aj_ can make the extrapolation 
of the A e ff(fc) plots unreliable. However, since the A e g (fc) data are really well described by a parabola, 
it might be true that the amplitude of such a correction term is small and the estimate Aj_ = 1.9723(90) 
is quite reasonable. In any case, this estimation shows a small deviation from the standard-theoretical 
picture, where 13. 1) should hold at small enough k with Aj_ = 2. This deviation can be indeed small at 
n = 10, since Aj_ — * 2 is expected in the limit n — ► oo, corresponding to the known behavior of the spheri- 
cal model 1181 . 

We have also attempted to evaluate the longitudinal exponent X\\ from the Gj (k) data within k> k* , 
where k* is indicated in figure[2]by a vertical dashed line. We have found that the longitudinal effective 
exponent, extracted from the data within [k,4k], can be perfectly approximated by a parabola. It leads to 
an estimate Ay = 0.85 ± 0.06. The error bars indicated here include a statistical standard error as well as 
a systematical error due to finite- ft effects. However, due to a rather large extrapolation gap (from « 1.17 
to ~ 0.85), we consider this estimation as a preliminary one. 
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Figure 4. The plots of the transverse effective exponent A e ff(fc), evaluated from the lnGx(k) vs Ink fits 
within [k,6k]. The results of the 0(10) model at ft = 3 are shown by solid circles, whereas those of the 
0(4) model at /5 = 1.1 — by empty circles. 

We have reexamined the largest-L (L = 350) and smallest-fa (ft = 0.0003125) data of the 0(4) model 
116] at p — 1.1 with an aim to evaluate the transverse effective exponent in the same way as for the 0(10) 
model. Like in the n — 10 case, we have verified that the thermodynamic limit at ft - +0 is practically 
reached in this estimation. Besides, we have found a surprising similarity of the A e ff(fc) plots, where the 
effective exponent in both cases was evaluated by fitting the G± (k) data within [k, 6k] (fits within [k, 4k] 
were used in 1 16]). As we can see from figure [4] the plot for the 0(4) model is systematically shifted 
down by an almost constant value relative to the same plot for the 0(10) model. This similarity might be 
partly caused by the fact that the values of spontaneous magnetization are rather similar in these two 
cases, i. e., M = M(+0) = 0.484475(48) for the 0(4) model at )3 = 1.1 and M a 0.467 (see section[5} for 
the 0(10) model at p = 3. The overall fit to 13.1) is less perfect for the 0(4) model as compared to the 
0(10) model. However, the systematical shift between two plots is very well approximated by a constant 
value within the statistical error bars for the eight smallest k data points in figurefj] It yields an estimate 
AAj_ = (Ax) n =io _ C^J-)n=4 - 0.0121(52), where (AjJ n =io and (Aj_) n=4 are the values of A± in n — 10 and 
n — 4 cases. According to the behavior of plots in figure |4]and this estimation, it is quite plausible that a 
transverse exponent Aj_ of the 0(10) model is somewhat larger than that of the 0(4) model. It is consistent 
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with the idea that 2- Aj_ decreases for large n and tends to zero at n — oo. This behavior is fully consistent 
with the predictions of |8], but not so well consistent with the standard theory, according to which Aj_ 
is always 2 and, therefore, AAj_ = is expected. According to our estimates (A_i_)„ =10 = 1.9723(90) and 
AAj_ = 0.0121(52), we have Aj_ = 1.960(10) for n = 4. It perfectly agrees with our earlier estimate Aj_ = 
1.955 ±0.020 111. 



4. The ratio universality test 

We have extended the MC analysis of our earlier data 11611 for the 0(4) model at two different cou- 
plings, p - 1.1 and p = 1.2, to test the expected (according to | 8]) universality of the ratio bM 2 /a 2 , dis- 
cussed already in the end of section [TJ According to <1.5) , dl.6t and <1.8> , the universality of bM 2 1 a 2 
implies that the quantity 

fc- rf M 2 Gii (k) 

tends to some universal constant at k — ► 0, i. e., limi?(k) = bM 2 la 2 . We have tested this property by 

fc-o 

comparing the R{k) plots at p = 1.1 and p = 1.2, where R{k) = R{\ k |) in the (100) direction. Note that the 
quantities in f4.lt are determined in the thermodynamic limit at fa — +0. We have verified that this limit is 
practically (within the statistical error bars) reached within k 3= A44 (where k( = 2n£/350) for the largest 
lattice size L = 350 and the smallest external fields h = 0.0003125 and h - 0.0004375 at which simulations 
were performed. The estimates of spontaneous magnetization obtained in (8J], i.e., M = 0.484475(48) at 
p = 1.1 and M = 0.560178(40) at jS = 1.2, are used here. The calculated plots are depicted in figured] The 
results for both h = 0.0003125 and h = 0.0004375 are available at p = 1.1. As we can see from figure[5] 
the corresponding two plots of R(k) (solid circles and diamonds) lie practically on top of each other, 
indicating that the finite-fa effects are negligibly small. The range h 3= 0.0004375 is considered for ft = 1.2 
in rial . Fortunately the finite-fa effects at f> — 1.2 are similar to those at p — 1.1, so that the estimate 
of R{k) at fa = 0.0004375 is valid at /3 = 1.2. The corresponding plot (empty circles) in figure [5] slightly 
deviates from the two plots at p — 1.1. However, all three plots merge within the statistical error bars at 
the smallest wave vector magnitudes k considered here. This confirms the expected universality of the 
ratio bM 2 1 a 2 . The plot of empty circles in figure [5] apparently saturates at a value about 0.166 for small 
wave vectors. Taking into account the two other plots, we can judge that 0.16 < R{0) < 0.18 most probably 
holds for the asymptotic value R(0) = lim i?(fc). Thus, we have an estimate R(0) = 0.17 ± 0.01. 




0.2 0.4 0.6 0.8 1 1.2 k 

Figure 5. The R(k) plots of the 3D 0(4) model, evaluated from the MC data for the lattice of size L = 350 at 
P = 1.1 and h = 0.0003125 (solid circles), j6 = 1.1 and h = 0.0004375 (diamonds), as well as at p = 1.2 and 
h = 0.0004375 (empty circles). 
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5. Spontaneous magnetization 

We estimated the spontaneous magnetization of the 3D O(10) model at p" = 3 based on our mag- 
netization data M{h,L) depending on h and L. We observed a rather fast convergence to the thermo- 
dynamic limit, e. g., M{h min ,L) = 0.31658(46), 0.45286(23), 0.470839(90), 0.471959(38), 0.472151(23), 
0.472148(16) at L = 32,64, 128, 192,256 and 350, respectively. According to this, we can take the largest- 
L value as a good approximation for M(h m i n ) = limi^ 00 M(/i m i n ,L). Using this method, we obtained 
M(Vin) = 0.472148(16), M{2h min ) = 0.4742786(98) and M(4/i min ) = 0.4772753(76). According to IPl , 
we fit these data to the ansatz M{h) - M(+G) + a\h p with p = 0.5211(69), evaluated from (1.9) by inserting 
here A± = 1.9723(90) obtained in section[3] It yields M = M(+0) = 0.467343(99). Assuming the standard- 
theoretical value p = 1/2, we obtain M = 0.467030(26). Fits to a refined ansatz M(h) = M(+0) + aih p + a 2 h 
yield M = 0.46711(14) at p = 0.5211(69) and M = 0.46696(14) at p = 1/2. Since we have only three data 
points for M(h), this can be considered as a raw estimation yielding M a 0.467. However, this estimation 
is accurate enough to see that p-3 corresponds to the ordered phase with M > 0. 

6. Conclusions 

In the actual work, the previous MC studies (l(J [3 of the transverse and longitudinal correlation 
functions in the 3D 0{n) models with n — 2 and n — 4 have been extended, including the n — 10 case 
(sections[2]and[3). It gives us an important information about the behavior of the exponent Aj. at large 
n. According to our MC analysis, a self-consistent (within the standard theory) estimation of A± for n = 
10 shows a small deviation from the standard-theoretical prediction Ax = 2, yielding Ax = 1.9723(90) 
(section[3}. The fact that this deviation is quite small can be well understood, since Ax — * 2 is expected at 
n — » oo according to the known results for the spherical model, corresponding to this limit. Comparing 
the plots of the effective transverse exponent at n — 10 and n — 4, it has been stated that these plots 
are surprisingly similar, i. e., only slightly shifted with respect to each other. The estimation of this shift 
suggests that the transverse exponent for n — 10 is larger than that for n — 4 by an amount of AAx = 
0.0121(52) (section [3). It is consistent with the idea that 2 - Ax decreases for large n and tends to zero 
at n — > oo. We have also verified and confirmed the expected universality of the ratio bM 2 1 a 2 for the 
0(4) model by analyzing the correlation functions at two different couplings, i. e., p 1 = 1.1 and p 1 = 1.2 
(section[4). 

The actual MC results are fully consistent with the predictions of the GFD theory |8] (see section [1} 
and not so well consistent with the standard theory, according to which Ax is always 2. 
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Ci/iHry/iflpHOCTi rojiflCToymBCbKHX MOfl b O(n) MOfle/iflx 

fl. Kayny3(P! P.B.H. Me/ibHii<P, fl. PiMcaHcP 

iHCTHTyr MaTeMaTHKH Ta KOMn'KiTepHwx HayK, YHiBepcnTeT /laTBii, 
6y/ibBap fl. PafiHica, 29, LV-1459 Pura, TlaTBm 

iHCTHTyT MaTeMaTUHHUx HayK Ta iHcfmpMaujCiHHX TexHO/iorifi, YHiBepcnTeT m. /lienan, 
By/i. /lie/ia, 14, LV-3401 f\ienan, /laTBia 

YHiBepcnTeT iM. Bi/ib<f>pefla /lopbe, BaTep/ioo, Omapio, KaHafla, N2L3C5 

y TpuBMMipHi/ix 0(n) MOfle/inx 3 n = 2,4,10 3flitiCHeHO aHa/ii3 MeTOflOM MoHTe Kap/io (MK) CHHry/wpHO- 
creti ro/iflCTOymBCbKi/ix MOfl fl/ia nonepenHOi i no3flOB>KHboT Kope/iflLjMHWx d>yHKi4iR, flKi noBOflHTb ce6e hk 
Gx(k) - ak~^ L i G||(k) = bk~^ y BnopaflKOBami/i cf>a3i npn k^O. HauiOK) MeTOK) e nepeBipnTi/i L^iKaBi Teo- 
peTMHHi nepefl6aHeHfl, 3riflHO hki/ix iHfleKCH X± i A|| e HeTpuBia/ibHWMH (3/2 < X± < 2 i < A| < 1 y Tpbox bh- 
Mipax) i KoecfiiujeHT bM 2 1 a 2 (fle M e cnoHTaHHOK) HaMamiHeHic™) e yHiBepca/ibHwtf. TpuBia/ibHi craHflapTHi 
TeopeTH4Hi 3Ha4eHHH e Xj_ = 2 i An = 1. Haiu nonepeflHifi MK aHa/ii3 flae X± = 1.955 ± 0.020 i A|| npn6/in3H0 
piBHe 0.9 fl/ifl 0(4) MOfle/ii. HeflaBHH MK oujHKa A|, flKa flonycKae nonpaBKn fl/ia CKefi/iiHra CTaHflapTHOi MOfle- 
n\, flae A|| = 0.69 + 0.10 ftfln 0(2) MOfle/ii. Tenep mu 3flMCHW/ii<i nofli6Hy MK oujHKy fl/ia 0(10) MOfle/ii, flKa flae 
Aj_ = 1.9723(90). Mu no6aHH/in, luo rpacpiK etfieKTHBHoro nonepenHoro iHfleKca fl/ifl 0(4) MOfle/ii e cucreMa- 
tmhho 3cyHy™Pi bhh3 no BiflHoweHHK) flo rpacfiiKa fl/w 0(10) MOfle/ii Ha AAx = 0.0121(52). I^e y3ro,q>KyeTbCfl 
3 flyMKOK), luo 2- Aj_ 3MeHi±iyeTbCfl pj\ft Be/inKnx n i npflMye flo Hy/ia npi/i n -» oo. Mm TaKOx nepeBipn/in i 
niflTBepfln^w oniKyBaHy ymBepca/ibHicrb bM 2 la 2 fl/ia 0(4) MOfle/ii, ppn aK01 6y/io 3flificHeHO cuMy/iflijii' npn 
flBOx pi3Hnx TeMnepaTypax. 

lOnowoBi c/iOBa: MOflejuoBaw-in Mome Kapno, n-KOMnoHGHTHi BeKiopni Moment, Kopenm4wn\ <pyHKU,i'i, 

CHHry/lfipHOCTi rOJlflCTOyHIBCbKI/IX MOfl 
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